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Abstract. In this paper, we define a new product over R°°, which allows us 
to obtain a group isomorphic to R* with the usual product. This operation 
unexpectedly offers an interpretation of the Redei rational functions, making more 
clear some of their properties, and leads to another product, which generates a 
group structure over the Pell hyperbola. Finally, we join together these results, 
in order to evaluate solutions of Pell equation in an original way. 



1. Introduction 

The Pell equation is 

x 2 — dy 2 = 1, 

where d is a given positive integer and x, y are unknown numbers, whose values we 
are seeking over the integers. If d is a perfect square, the Pell equation has only the 
trivial solution (1,0). The interesting case is when d is not a square, because then 
the Pell equation has infinite integer solutions. For example, knowing the minimal 
non trivial solution (xi,yi) with x\ and y\ positive integers, we can evaluate 

[x\ + yiVd) n = x n + y n Vd, Vn > 

where (x n ,y n ) is also a solution of the Pell equation (cf. [5]). The Pell equation has 
ancient origins (for further references, see pQ, [14]) dating back to Archimedes Cattle 
Problem, but it is not known if Archimedes was able to solve it. Surely the first 
mathematician who found a method for solving the Pell equation was the Indian 
Brahmagupta. If you know two solutions (xi,y±) and (#2,2/2) of the Pell equation, 
with the Brahmaguptha method you find that 

(xix 2 + dyiy 2 ,xiy 2 + yix 2 ), (x x x 2 - dy x y 2 , x x y 2 + y x x 2 ) 

are also solutions. Nowadays the most common method of resolution involves con- 
tinued fractions. A continued fraction is a representation of a real number a through 
a sequence of integers as follows: 

1 

a = clq H 



ai + 



1 



1 



a 2 + 

a 3 H 

where the integers ao,ai, ••• can be evaluated with the recurrence relations 
Ofc = [ak] 

1 fc = 0,1,2, ... 

ak+i = if ak is not an integer 

ak~ a k 

1 
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for ao = a (cf. A continued fraction can be expressed in a compact way using 

the notation [ao, ai, 02, 03, ...]. The finite continued fraction 

[ao,...,a n ] = — , rc = 0, 1,2, ... 

Qn 

is a rational number and is called the n-th convergent of [ao, ai, a>2, o>3, ■■■]■ The 
continued fraction expansion of Vd is periodic and it has the form (cf. |10j ) 

Vd= [a ,ai, ...,a L _i,2a ] . 

It is possible to prove that the minimal solution of the Pell equation is 

= (PL-i,qL-i), L even 

(xi,yi) = (p2£-l,02£-i), odd, 

where 

PL-l 



QL-1 
V2L-1 



[a ,...,a L -i] 
[ao, a2L-i] • 

Q2L-1 

Furthermore all solutions of the Pell equation are of the form (p n ^_i, q n i-i), when 
L is even, and (p2nL-i, Q2nL-i) when L is odd, for n = 1, 2, ...( cf. [10]). 
The geometrical locus containing all solutions of the Pell equation is the Pell hyper- 
bola 

H d = {(x,y) £R 2 :x 2 -dy 2 = 1} . 

The fascinating aspect of working with the Pell hyperbola is a group law that we 
can give over H^. In fact, for every pair of points P,Q £ H^, we can define their 
product as the intersection between and a line through (1,0) parallel to the line 
PQ. This geometric construction is essentially the same presented in the book of 
Veblen (cf. [13]) as a general product between points on a line ([13] p. 144), which 
can be defined in a similar way over conies ([13], p. 231-232), in order to obtain 
a group structure. Algebraically (see, e.g., [5]) this product between two points 
(xi,yi) and (#2,2/2) 1S given by the point 

{xxx-2, + dyiy 2l y\x 2 + xxy 2 ) ■ 

This geometric group law over the Pell hyperbola, similar to that one over an el- 
liptic curve, has the advantage of a nice algebraic expression, which yields greatly 
simplified formulas (cf. [6]). In this paper we will see how it is possible to obtain 
this product in an original way, starting essentially from a transform over RU {00}. 
Then we will use the Redei rational functions (cf. [11]) in order to evaluate the 
powers of points over Hj. 

The Redei rational functions arise from the development of (z + vd) n , where z is 
an integer and where d is a nonsquare positive integer. One can write 

(z + Vd) n = N n (d,z) + D n (d,z)Vd , (1.1) 

where 

[n/2] [n/2] 
N n (d,z) = £ (;>^- 2fc , D n (d,z) = £ L n W-»-i. 
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The Redei rational functions Q n (d,z) are defined by 

Qn(d,z) = ^ n{d,Z) Vn>l. (1.2) 
D n (d,z) 

It is well-known their multiplicative property 

Q nm (d,z) = Q n (d,Q m (d,z)) , 

for any couple of indexes n, m. Thus the Redei functions are closed with respect to 
composition and satisfy the commutative property 

Q n (d,Q m (d,z)) = Q m (d,Q n {d,z)) . 

The multiplicative property of the Redei functions has several applications. For 
example, this property has been exploited in order to create a public key crypto- 
graphic system (cf. [9]). The Redei rational functions reveal their utility in several 
other fields. Given a finite field ¥ q , of order q, and Vd ¥ q , then Q n (d,z) is a 
permutation of ¥ q if and only if (n,q + 1) = 1 (see [7j, P- 44). Another recent 
application of these functions is in finding a new bound for multiplicative character 
sums of nonlinear recurring sequences (cf. [1]). Moreover, they can be used in order 
to generate pseudorandom sequences (cf. [12]). In this paper we will see a totally 
different approach for the Redei rational functions, studying an original connection 
between them and continued fractions. In particular we will see how we can apply 
Redei rational functions in order to generate solutions of the Pell equation. This 
research starts from the fact that the sequence of functions Q n (d,z), by definition, 
converges to y/d, i.e., Redei functions are rational approximations of Vd, for any pa- 
rameter z. We know that the best rational approximations of an irrational number 
are provided by the convergents of its continued fraction (cf. [10J). So it would be 
beautiful to find a parameter z in order to obtain some convergents from Q n (d, z). 
In this paper we find the parameter z such that Q2n(d,z) correspond to all the 
convergents of the continued fraction of Vd, leading to the solutions of the Pell 
equation. This method uses the Redei rational functions in a totally different field 
from the classic ones and, considering the fast evaluation of Redei rational functions 
(cf. [8]), it allows to generate solutions of the Pell equation in a rapid way. 

2. A NEW OPERATION OVER R°° 

We consider a transform, defined over R°° = M. U {oo} which induces a natural 
product over We show how this operation is strictly related to the Redei 

rational functions (cf. [11]). revealing a new interesting point of view to understand 
their multiplicative property. 

Definition 2.1. For any positive real number d, we define the transform 

Pd : R°° R°° 

Pd {x) = ^—Vd . (2.1) 
x — 1 



As immediate consequences of (|2.ip . we have that 

p d (l) = oo, Pd {0) = -Vd, p d (oo) = Vd . (2.2) 
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Moreover, we can easily find the inverse of pd- 

Pd \x) = (2.3) 

x — v a 

A more important achievement is the product 0^, induced by pd over 



P X. 



i i d + xv 

xQ d y = Pd {p- d l (x) Pd l (y)) = —— , (2.4) 

x + y 

which is surely associative and commutative. Furthermore, comparing with the 
usual product over R, \fd plays the role of oo and —yd the role of 0, since we have 
the following relations for all x S R°° 

V~dQdX = — ~t- = Vd and — VdQdX = —Vd. 
V d + x 



These properties, together with (|2.2p . suggest that this new product induces a group 
structure on the set 

P d := R°° - {±Vd} . 

Over (P d ,Qd), o° is the identity with respect to and any element x has the 
unique inverse — x since 

x Q d oo = x, x Qd (—x) = oo, Vx G P d ■ 

We immediately have the following. 

Theorem 2.2. (Pd,Q d ) is a commutative group. The transform pd is an isomor- 
phism from the group (R*, •), of nonzero real numbers under ordinary multiplication, 
into the group (Pd,Qd)- 

Remark 2.3. If we consider two positive real numbers d, e, we also have an imme- 
diate isomorphism between (Pd,Qd) an d (P e ,& e ) : 

<p:P e ^Pd 

I* 

(p '. x i — y X\ — . 
V e 

We now show a wonderful relation relating the product (|2.4p and the Redei ra- 
tional functions (cf. It is well-known that for any pair of indices m and 

n, 

Qnm(d, z) = Q n (d, Q m (d, z)) (2.5) 

but what can we say about Q n+m (d, z)l Using product fj2.4|) the answer is simple 
and furthermore we can use this fact in order to obtain the multiplicative property 
in a more explicative way than the usual one (see [7], p. 22-23). 

Proposition 2.4. With the notation introduced above 

Qn+m (d, z) = Q n (d, z) Q d Qm(d, z). 

Proof. For the sake of simplicity, here we omit the variables z and d, so 

N n (d, z) = N n , D n (d, z) = D n , Q n (d, z) = Q n . 
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Using a matricial approach (cf. [3]) we have 

z d\ n _(N n dD r 

1 z ~\D n N n 



(2.6) 



so 

\ ( N n dD n \ ( N m dD m 
, D n+m N n+m ) ~ \D n N n ) \D m N m 

Comparing the resulting matrix on the right with the one on the left, we finally 
obtain the relations 

(N n+m = N n N m + dD n D m ^ 

[Ai+m = D n N m + N n D m . 

Now the proof is straightforward using (|2.7j) and definition of Redei rational func- 
tions (OP : 

N n N„ 



d + QnQm d +Tt-T^ N n N m + dD n D Ti 

N n I N„ 

D n ^ Dr, 



VnW^m- fl , - Nn Nm ~ AT f) i AT tj ~ ^n+vn 



□ 



Remark 2.5. It is interesting to consider the additive property of the Redei ra- 
tional functions because it clarifies the meaning of the multiplicative property and 
makes this relation simpler to use. Furthermore, we can relate these functions to 
the integers with the operations of sum and product, extending their definition to 
negative indexes n through this matricial approach. Indeed the set of the Redei ra- 
tional functions {Q n (d, z),\/n G Z} with the operations Q d and o (the composition 
of functions) is a ring isomorphic to (Z, +, •). 

The previous proposition is a powerful tool, which enables us to evaluate the 
powers of elements with respect to the product (12. 4p . 



Corollary 2.6. Let z nQ d = z Q d • • • 0^ z be the n-th power of z with respect to the 

V v ' 

n 

product \2.J$. Then 



z n ^=Q n (d,z) . 
Proof. From (jl.ip and fjl .2[) . when n = 1 

z = Qi(d,z) . 

Therefore 

z n Q d = zQ d ---Q d z = Q 1 ( y d,z)Q d ---Q d Q ± (d, z) = Q 1+ ... +1 (d, z) = Q n (d, z). 

□ 

This Corollary shines a light on the multiplicative property (|2.5p . which holds 
because Q n (d,z) is essentially a power of an element with respect to the product 
H2T 



Q n (d,Q m (d,z)) = (Q m (d,z)) n ° d = (z m °i) n °i = z nm °* = Qnm{d,z). 

This proof for the multiplicative property allows us to consider the Redei rational 
functions as a power of an element and in this way they are easier to use, as we will 
see later. 
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3. A COMMUTATIVE GROUP OVER PELL HYPERBOLA 

The group (Pd,Qd), related to a nonsquare positive integer d, reveals some im- 
portant connections to the Pell hyperbola 

H d = {(x,y) £R 2 /x 2 -dy 2 = 1}, 

with asymptotes 

x 

y = ± 7=d- 

The Pell hyperbola Hd can be parametrically represented using the line 

y = -(x + l). (3.1) 

m 

Here we use, as parameter to, the cotangent of the angle formed with the x-axis, 
instead of the usual tangent. We choose the cotangent as parameter because it is the 
only choice which allows us to have a correspondence between Pd and Hd- Indeed 
in this case oo (the identity of Pd) will correspond to (1,0) (the identity of Hd). 
Moreover, the two points at the infinity are obtained by the parameters ±vd. 
Now, for all m ^ iy/d, the line (I3.ip intersects Hd in only two points: ( — 1,0) and 
another one. Thus it seems natural to visualize Pd as the parametric line. Moreover, 
the bijection 

£d : Pd — > H A 
' m 2 + d 2m 



v to z — d m L — d 



yields a solution to the system 



x 2 - dy 2 = 1 



1 , . (3-2) 
— (x + 1) . 



m 

In fact, we find two solutions (—1, 0), , 1 ^I n i d J ■ The bijection ed maps into 

(—1,0), oo into (1,0), and, using (|3. If) . the inverse of td can be seen to be 

t : H d -> P d 

. l + x 
t : (x,y) i y . 

y 

We observe that the map r has no explicit dependence on d, but this dependence is 
implicit when we consider a point (x, y) E H d - Since we have a bijection from Pd to 
Hd, we can transform the product over P d into a product 0# over H d as 

(s, *) ®H (u, v) = e d (r(s, t) @ d t(u, v)), V(s, t), (n, v) G H d . 

Using the definitions of e and r: 

/l + s\ f 1 + u\ l + s + u + su + dtv 
r(s,t) Q d T(u,v) = -—) Q d 



t ) \ v J t + tu + V + sv 

and, using the relations s 2 — dt 2 = 1 and u 2 — dv 2 = 1, finally we have 

(s, t) Qh (u, v) = (su + dtv, tu + sv). (3-3) 
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The operation constructed on Hd has all the good properties of 0^. Therefore we 
have the following 

Theorem 3.1. {H^Qh) is isomorphic to (Pd,Qd) and (Hd,Qi{) is a commutative 
group. 

As we have seen in the introduction, the product f)3.3[) is a classical product used 
to construct a group over a conic, but here we have seen how it is possible to obtain 
this product in an original way, starting essentially from a transform over R°°. 

Remark 3.2. Hd is isomorphic to the group of unit norm elements in Q(\/d)- 
Remembering that 

Q(Vd) = {s + tVd: s,t G Q}, 
the natural product between two elements of Q(^/d) is 

(s + tVd)(u + vVd) = (su + dtv) + {tu + sv)\Q. 

Remark 3.3. We can easily observe that the products Qd and Oh> an d the transfor- 
mations ed, t, all involve only rational operations. If we consider P^ = Q U {oo} = 
Q°° andWd = {(x,y) 6 Hd/x,y G Q}, then we retrieve again €d,r, as isomorphisms 
between and Md- Furthermore does not depend on d and (P^, Qd) — (Q*i ■)■ 

4. Generating solutions of the Pell equation via Redei rational 

functions 

We are ready to put together all the results obtained in the previous sections. The 
main purpose of this work is to reveal a beautiful connection among Redei rational 
functions, the product Oh and solutions of the Pell equation. 

The matrix representation of N n (d, z), D n (d, z), introduced in Proposition 12.41 al- 
lows us to consider 

(N n (d,z))+™=W(l,z,2z,z 2 -d) 

{D n {d,z))+™=W{^l,2z,z 2 -d) . 

Here (a n )^^ = W(a, b, h, k) indicates the linear recurrent sequence of order 2, with 
initial conditions a,b and characteristic polynomial t 2 — ht + k, i.e., 

ao = a 

< a\ = b 

a n = ha n -\ — ka n -2 Vn > 2 . 

Now for any point (x, y) £ Hd we set 

(x n ,y n ) = {x,y) n °n 7 

where 

0o,yo) = (1,0) and (x 1 ,y 1 ) = (x,y). 

We know that 

(x + yVd) n = x n + y n Vd , 
and so (x n ), (y n ), as sequences, recur with polynomial t 2 — 2xt + 1. We can easily 
observe that (x,y) n@ H = (F n (x),yG n (x)) and 

(F n (x))+™=W(l,x,2x,l) 

(G n (x))+r = W(0,l,2x,l) . 
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Now, comparing the recurrence of F n (x) and N n (d, z), we can see that they coincide 
when 



introduced by L. E. Dickson over finite fields. He studied when they give permu- 
tation of the elements of the finite fields. Furthermore, Dickson polynomials are 
related to the classical Chebyshev polynomials and they are used in several areas, 
like cryptography, pseudoprimality testing, in order to construct irreducible poly- 
nomials over finite fields and in many other applications of the number theory. The 
Redei rational functions are related to the Dickson polynomials: 



When a = 1, the sequence (g n (l,x) = g n ( x ))n=o recurs with characteristic polyno- 
mial t 2 — xt + 1 and we can observe that 



Since (x,y) nm@ H = ((x, y) n ®H )"% ; F n (F m (x)) = F nm (x) and we have retrieved 
the multiplicative property of g n (x). Furthermore we have the explicit formula for 
F n (x): 



z = x and z 2 — d = 1 , 
i.e., d = z 2 — 1 = x 2 — 1. We immediately have the following 

Proposition 4.1. With notation introduced above 

F n (x) = N n (x 2 — 1 , x) and G n (x) = D n (x 2 — 1, x). 
Let us recall Dickson polynomials (cf. [2]): 




2N n (d,z) = g n (2z,z 2 -d). 





Q n {d, a) Q d Q n {d, b). 



We have 




and when a = b = z 




d + Q n (d,z) 2 



2Q n (d,z) 



Q n (d, z) Q d Q n (d, z) = Q 2n (d, z). 



Given (x, y) S H^, we obtain 
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Since (r(a;,y)) n0 d = T ((x, y) n ©» ) , 

b " ( ' V) r{{F n {x),yG n {x))) = Q n (d 1 + ' 



and thus 

'l + x\ n0d 1 + A" n (x 2 - 



Proposition 4.2. For any (rr,y) G Hj, 

2n 



J yD n (x 2 -l,x) ' 



l + xV n0d F n (x) A?" n (x 2 -l,x) 



y J yG n (x) yD n (x 2 -l,x)' 

Proof. Here, we write N n ,D n , instead of N n (x 2 — 1, x), D n (x 2 — l,x). We want to 
prove the equality 

1 + F 2n (x) F n (x) 



yG2n(x) yG n (x) ' 

and so we consider 

1 + N 2n N n _ D n + D n N 2n - N n D 2n 
yD 2n yD n yD 2n D n 
As a consequence of (|2.7f) we have 

N 2n {d, z) = N*(d, z) + dD 2 n {d, z) 
D 2n (d,z) =2D n (d,z)N n (d,z) , 

and 

l + iV 2n N n D n {l - N 2 + (x 2 - l)D 2 n ) 



yD 2n yD n yD 2n D n 

In fact, using (|2.6[) 

N 2 (d,z)-D 2 n (d,z) = (z 2 -d) n , 

and in our case 

N 2 (x 2 - l,x) - (x 2 - l)D 2 n (x 2 -l,x) = 1. 



0. 



□ 



Corollary 4.3. With the notation at the beginning of Section 3, for any {x\, y\) £ 
Hd, i.e., such that x\ — dy\ = 1, we have 

( X\ + 1\ x n 
Q2n d, = — . 

V yi / Vn 

We want to apply all these tools, and in particular the Redei rational functions, 
to the solutions of the Pell equation x 2 — dy 2 = 1. In what follows, we will use Redei 
rational functions to find convergents of continued fractions which provide solutions 
of the Pell equation. 

Pn 

Let (x\,yi) be the minimal integer solution of the Pell equation and let — , for 

q n 

n = 0,1, 2, be the convergents of the continued fraction of ^/d and L the length 
of the period. As we have seen in the introduction, it is well-known that {x\,y\j = 
(pL-i,qL-i) when L is even and (x\,yi) = (p 2 L-i, <?2L-i) when L is odd (see [TP]). 
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Furthermore, all the solutions of the Pell equation are (p n L-i, q n L-i), when L is 
even, and (p2nL-i,Q2nL-i), when L is odd, for n = 1, 2, ... . Then we have 

(xi,yi) n0 « = (x n ,y n ) = (pnL-i, qnL-i) 

or 

(xi,yi) nQH = (x n ,y n ) = (P2nL-l,q2nL-l) , 

when L is even or odd respectively, and by the last Corollary: 

Q2n a, = , L even 

V VI J QnL-1 



Q2n (d, 1 + ) = P2nL - , L odd . 

V 2/1 / <72nL-l 

In other words, using Redei rational functions, we can evaluate all the convergents 
of the continued fraction of \fd leading to the solutions of the Pell equation. 
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